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Abstract 

We apply the methodology of optimal control theory to the problem of implementing quantum 
gates in continuous variable systems with quadratic Hamiltonians. We demonstrate that it is 
possible to define a fidelity measure for continuous variable (CV) gate optimization that is 
devoid of traps, such that the search for optimal control fields using local algorithms will not 
be hindered. The optimal control of several quantum computing gates, as well as that of 
algorithms composed of these primitives, is investigated using several typical physical models 
and compared for discrete and continuous quantum systems. Numerical simulations indicate 
that the optimization of generic CV quantum gates is inherently more expensive than that of 
generic discrete variable quantum gates, and that the exact-time controllability of CV systems 
plays an important role in determining the maximum achievable gate fidelity. The resulting 
optimal control fields typically display more complicated Fourier spectra that suggest a richer 
variety of possible control mechanisms. Moreover, the ability to control interactions between 
qunits is important for delimiting the total control fluence. The comparative ability of current 
experimental protocols to implement such time-dependent controls may help determine which 
physical incarnations of CV quantum information processing will be the easiest to implement 
with optimal fidelity. 
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I. INTRODUCTION 



The optimal control of quantum dynamical systems has become a subject of intense 
interest in chemistry, physics and most recently, information theory jl, 2]. Over the past 
several years, it has become clear that the physical implementation of logical gates in 
quantum information processing (QIP) may be facilitated by using the methods of optimal 
control theory (OCT) [si, 0, [5, 6, 7, ^j]. When implementing a quantum logic gate through 
OCT, the distance between the real and ideal quantum unitary transformation can be 

n 

used as a metric for assessing optimality [3J. The map between admissible controls and 
associated values of this distance is referred to as the control landscape; extremal control 
solutions to the problem correspond to critical points of this map. These landscapes were 
recently shown to universally possess very simple critical topologies for finite-dimensional 
quantum gates, with no suboptimal traps impe ding optimal searches, irrespective of the 



lot 



system Hamiltonian for controllable systems 

Prior work on the implementation of quantum gates using OCT has been directed to- 
ward discrete QIP, as hypothetically carried out on the so-called quantum spin computers 
originally discussed by Benioff and Feynman 
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12j. These computers, in which infor- 



mation is carried as quantum bits (qubits) |l3|] encoded in discrete systems like electron 
spins or two-level atoms, are the quantum version of digital classical computers. Classical 
information can be carried by either a discrete (digital) signal or a continuous (analog) 
signal. Quantum information can also be carried by continuous (infinite-dimensional) 
systems, such as a harmonic oscillator, rotor, or modes of the electromagnetic field [l4|. 
Quantum information processing over continuous variables (CV) can be thought of as the 
systematic creation and manipulation of continuous quantum bits, or qunits 14]. 

Importantly, CV QIP may be less susceptible to drift than its classical counterpart. 
Cleverly encoded quantum states can be restandardized and protected from the grad- 
ual accumulation of small errors, or from the destructive effects of decoherence fis| . 
Moreover, compared to discrete QIP, continuous-variable QIP has several practical ad- 
vantages, originating for example in the high bandwidth of continuous degrees of freedom, 
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that have spurred substantial interest in its generic properties. Significant advances have 
recently been made in the experimental implementation of continuous QIP, including the 
demonstration of quantum teleportation over continuous variables 16|. As experimen- 
tal methodologies improve, it becomes important to consider how such implementations 
could be enhanced through the systematic application of OCT. 

Early control studies on continuous quantum systems focused on the manipulation of 
quantum scattering states in bond-selective control (i.e., dissociation or association of 
atoms) of molecular systems [ill, Q. Such systems were shown to be associated with 
dynamical symmetries represented by noncompact Lie groups with infinite-dimensional 
unitary representations [19|. A criterion for approximate strong controllability was given, 
showing that such systems, which possess an uncountable number of levels, can be well 
manipulated using a finite number of control fields. 

Here, we examine OCT problems pertaining to an important class of continuous quan- 
tum gates, namely those that can be represented as symplectic transformations of the 
quadrature vectors in the Heisenberg picture. This gate set is referred to as the Clifford 



group 



20] . Although universal quantum computation over continuous variables requires 



higher-order nonlinear operational gates, networks using only Clifford group gates have 
numerous important applications in the area of quantum communication, and in fact, 
these gates are in many ways easier to implement over continuous variables than over 
discrete variables 21[. For example, quantum error correcting codes, which are essential 



for overcoming the effects of errors and decoherence, use only Clifford group gates for en- 
coding and decoding. Other important protocols in QIP, like quantum teleportation, also 
rely solely on Clifford group gates and related measurements. The Clifford group gates 
are sufficient to represent any CV quantum computation that can be efficiently simulated 
on a classical computer. As such, the symplectic gate formalism also has applications in 
the context of reversible analog classical computation. 

In this paper, we investigate the implementation of CV quantum gates via OCT. We 
carry out this analysis in two stages. First, we analyze the topology of the optimal control 
landscape for symplectic gate fidelity upon the assumption of full controllability of the 
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underlying control systems, which is demonstrated to be free of local traps that might 
otherwise impede the optimization process. This analysis is Hamiltonian-independent 
except for the assumption of full controllability and therefore the conclusions are generic 
for CV quantum computation. Next, we carry out numerical OCT calculations for various 
CV gates, comparing to the analogous discrete variable gates with two physical models. 
We identify characteristic differences between these two problems both in terms of control 
optimization efficiency and the complexity of the associated optimal Hamiltonians. 

The paper is organized as follows. Section II summarizes the symplectic geometry 
arising in quantum linear optics and provides preliminaries on the symplectic gates ap- 
plied in CV QIP. Section III analyzes the landscape topology for symplectic gate fidelity 
and its restriction on the maximal compact subgroup of the symplectic group. In Sec- 
tion IV, we discuss several physical models for CV QIP and compare their properties as 
control systems. In Section V, we carry out OCT calculations for specific CV gates and 
algorithms using these model Hamiltonians, comparing with the corresponding problems 
for discrete quantum gates. Finally, Section VI draws general conclusions regarding OCT 
for CV quantum gates versus discrete quantum gates. 



II. SYMPLECTIC GEOMETRY AND SYMPLECTIC GATES 



Formulation of optimal control theory for CV gates is simplified by framing the time 



evolution of the canonical observable operators in terms of group theory [22|. For example, 
suppose that the system is be realized as a quantized electromagnetic field. Let fij and a\ 
represent the creation and annihilation operators corresponding to the i-th mode of the 
field. These operators are related to the position and momentum operators by 

ft = + p l = ^(a, - fit) 

Then the time evolution of the vector of quadratures is associated with a time 
evolution of the 2iV-dimensional vector of creation and annihilation operators £ = 
(fiij ■ ■ ■ > %, • • • , ajy). 
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The system Hamiltonian H(t) generates a one-parameter family of evolution transfor- 
mations U(t) on the Hilbert space Ti that obeys the Schrodinger equation 

where the parameter is the time, and H is assumed to be a quadratic Hamilto- 
nian that takes on various forms depending on the nature of the coupling between 
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the oscillator modes [23(. Denote the vector of quadratures of quantum observables 
z = (<Zi> ■ ■ • , QN'iVii ' ' ' iPn) T ■ The evolution propagator transforms the quadrature vec- 
tors linearly through 

U: z a -> U^(t)z a U(t) = J2 S ^(t)z p , 

where the matrix S(t) is an element of the symplectic group Sp(2iV, R). The symplectic 
group is defined as the set of 2N x 2N matrices that satisfies S T JS = J, where 

In 

-In 

Thus, the matrix S captures the Heisenberg equations of motion for the operators Zi, and 
the unitary propagator U forms the metaplectic unitary representation of S in Sp(2iV, M). 
Correspondingly, let H — {hij} — H T be the matrix representation of the Hamiltonian, 
i.e., H(t) = YttijhjiteWji which belongs to the Lie algebra sp(2JV,R) of Sp(2iV,IR) = 
{JH | H T = H} 2J| with dimension being N(2N + 1). Through this representation, 
the symplectic matrix S(t) associated with U(t) follows a classical Hamiltonian evolution 
equation 

^ = JH(t)S(t). (2) 

The infinite dimensional unitary operator U(t) in ([1]) carries the (metaplectic) unitary 
representation of these symplectic transformations on the Hilbert space 7i of quantum 



states in the Schrodinger picture [23] 



Another important class of transformations are the displacements that shift q and p 
by constants: 

q — ► q + a, p — > p + b. 



The combination of displacements with homogeneous symplectic transformations forms 
the inhomogeneous or affine symplectic group ISp(2iV, R) in the form of 

S c = ( S ^ C \ S E Sp(2iV, R), c G M 2N , 
( z 

which acts on the extended phase vector 

V 1 

Symplectic operations executed by CV quantum computers are particularly important 
because they correspond to information processing tasks for which these computers are 
expected to outperform their discrete counterparts. In the context of quantum optics, 
which is the basis of most proposed schemes for CV quantum computation, they re- 
quire only linear optics and squeezing; as such, these gates may be fairly straightforward 



to implement [21]. Theoretically, the generalized Gottesman-Knill (GK) theorem pro- 
vides a valuable tool for assessing the classical complexity of a given continuous quantum 
information process [2^] . It states that any quantum algorithm that initiates in the com- 
putational basis and employs only the restricted class of affine symplectic gates, along 
with projective measurements in the computational basis, can be efficiently simulated 
on a classical computer. This computational model is often described within the stabi- 
lizer formalism 
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20| as Clifford group computation. To achieve universal continuous 
quantum computation, it is necessary to introduce additional operations (corresponding 
to elements of the nonlinear symplectic diffeomorphism group), afforded in the quantum 
optics laboratory by the ability to count photons. 

Quantum optical control components can include linear elements (such as beam split- 
ters, mirrors, and half- wave plates), nonlinear elements (such as squeezers, parametric 
amplifiers and down-converters) or a combination thereof. The symplectic transforma- 
tions that can be performed by linear optics consist of the inhomogeneous displacement 
transformation and the maximal compact subgroup OSp(2iV, R) of orthogonal symplectic 
matrices, which preserve the total photon number 

JV 



n 

i=i 
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This subgroup is isomorphic to the unitary group U (N) via the correspondence: 

X-iYeU(N) - S= ^j. (3) 

Squeezing operators (also referred to as active transformations) cannot be implemented 
using linear optics; they fall into the noncompact portion of Sp(2iV, R), and correspond to 
photon non-conserving transformations. These operators rescale canonical operators by 
a real number A along one axis in the quadrature plane, and by A -1 along the conjugate 
axis: 

q — > Xq, p — > X^p. 

There exists a set of universal symplectic gates whose combinations may realize arbi- 
trary Clifford gates in ISp(2iV, R). A well-known choice of the universal gate set consists 
of the Pauli operators, the Fourier gate, phase gate, and the SUM gate. The Pauli 
operators perform the phase displacements 

X(qi) = exp(ihqipi), Zfa) = exp^ftp^), 

whose symplectic representations are 





(l q\ 




( 1 \ 






X(q) = 


10 


, Z(p) = 


Ol-P 


G ISp(2iV,R). 


(4) 
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The one-qunit Fourier transform is the CV analog of the discrete Hadamard gate: 

(j)-(i 

This action can be represented by a 3 x 3 affine symplectic matrix 

/ 1 0\ 

F= -1 G ISp(2iV,R). (5) 
V o \) 

The phase gate is the analog of the discrete variable phase gate 



P(rj) = exp ( — r]q 



pi \ p — rjq 



Unlike the other gates, P is a function of a real parameter, and can be represented by 
the matrix 



(6) 



/ 1 o\ 

P(V) = -V 1 G ISp(2iV,R). 
V o \) 

Finally, the SUM gate acts on a two-qunit system where qunit 1 is said to be the 
control and qunit 2 is said to be the target [20j, [25(, and it carries out the following 



transformations on the canonical observable operators: 

SUM : q x -»• q 2 -> q x + q 2 , Pi -»• Pi - p 2 , P2 



P2- 



This is the continuous-variable analog of the discrete CNOT gate and its unitary repre- 
sentation is 



SUM = exp | -<i. 
Therefore, the associated symplectic representation is 



SUM 



/ 1 





0\ 


1 1 








1 


-1 








1 





^0 








e ISp(2iV,E). 



(7) 



III. CONTROL LANDSCAPE TOPOLOGY FOR SYMPLECTIC GATE FI- 
DELITY 



In this work, we are concerned with the optimal control problem of identifying the 
time-dependent functional form of the Hamiltonian that maximizes the fidelity of a sym- 
plectic (CV) gate at a fixed time tf, with a particular focus on the convergence of search 
algorithms for this problem. Prior work 26) has begun to examine the question of con- 
structing minimal time quantum circuits for a given symplectic gate from restricted con- 
trol Hamiltonians. By contrast, we are interested in the general problem of gate control 
for arbitrary quadratic Hamiltonians and arbitrary final times. This problem must be 
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solved through computational search over the space of time-dependent quadratic coupling 
Hamiltonians that minimize the distance to the target gate. In continuous quantum com- 
putation, we are primarily interested in the set of gates that can be efficiently simulated 
classically, i.e., those that form the inhomogeneous symplectic group ISp(2iV, R). Since 
the (quantum) symplectic gate U is a faithful unitary representation of a symplectic ma- 
trix S, it is reasonable and convenient to define the gate fidelity analogously to that for 
discrete gates as 

J[C(-),t f ] =Tr(S-W) T (S-W) + (s-w) T (s-w), S s G ISp(2iV, R), (8) 

where W w (including a homogeneous transformation W and a u>-translation) is the finite- 
dimensional representation of the target quantum propagator, and S s is the representation 
of the system propagator U(t) as an implicit function of C(-). Since the control of the 
inhomogeneous part via linear couplings has been extensively studied in the literature 



271 . |28I ]. we are primarily concerned with the control of the homogeneous part S and 
will always omit the inhomogeneous part (in practice, by switching off the corresponding 
linear interactions). 

The unitary propagator is an implicit function of the control field C(t) through the 
controlled Schrodinger equation 

dU(t) i 



dt h 



U(t), (9) 



where H is the (quadratic) internal Hamiltonian and Hi is the interaction Hamiltonian 
steered by a control field Cj(t) to couple the internal degrees of freedoms. Hence, the 
optimization problem is formally defined on the space of time-dependent control fields 
subject to the dynamical constraint of the Schrodinger equation (jHJ), which is equivalent 
to the following dynamical constraint on S: 



dS(t) _ j 



dt 



S(t). (10) 



H (t) + J2Q(t)H 

i=l 

Within the framework of quantum optics, a compact symplectic gate corresponds to a 
control implemented via linear optics, and a noncompact gate corresponds to a control 
involving implementations of squeezing. 



As described in prior work [9|, [29| , any search for the optimal control fields minimiz- 
ing such a cost functional for a given physical system will traverse a so-called control 
landscape, defined as the map between admissible control fields {Ci(t)} and the associ- 
ated values of J . In order to facilitate our comparison of the optimization efficiencies 
of discrete and CV quantum gates, it is useful to acquire the critical topology for the 
latter problem, i.e., the distribution of all possible critical points, including suboptima 
that may serve as undesirable attractors for the search trajectory. Since the cost func- 
tion is a complete function of the propagator S, it would be convenient to investigate 
the landscape topology on Sp(2iV, R). According to 0, 31|, the optimality status (i.e., 
minimum, maximum or saddle point) of a critical control field {Cj(i)} is equivalent to 
that of the resulting gate S on Sp(2iV, R), provided that the map from {Ci(t)} to S is 
locally surjective at {Cj(i)} on the set of realizable gates at the time tf. Such controls are 
conventionally called regular extremal controls. Vanishing of the gradient VJ7"[C(-)] may 
also be caused by singularities of the map from {Cj(i)} to S. These critical points, which 
are independent of the choice of the cost function, are called singular extremal controls 



32J. For simplicity, only regular extremals are considered, and the singular extremals 



will be studied in the future. Actually, in most cases of this paper the local optima are 



not singular [32]. 

In addition, we assume that the system is fully controllable at the final time tf, i.e., 
any Clifford group element can be achieved by some cleverly designed control functions 
{Cj(t)}, so that the system is capable of achieving arbitrary Clifford group computations. 
The critical topology problem can then be reduced from the (infinite-dimensional) do- 
main of control fields onto the symplectic group. Generally, a fundamental requirement 
for the controllability of systems evolving on Lie groups!^ is the rank condition, i.e., the 
condition that the Lie algebra spanned by H Q , Hi, ■ ■ ■ , H m , and their commutators such 
as [H , H], [H , [H , Hi}], [Hi, [H , Hj]], etc., equals the Lie algebra of the Lie group. This 
condition has been proven to be sufficient for the unitary propagators of discrete-state 



quantum systems 



and for these systems there exists a positive time Tq such that 



arbitrary gate can be achieved exactly at any positive time larger than Tq. For systems 
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on noncompact symplectic groups, the rank condition is only sufficient when H is com- 
pact, and there is no guarantee of exact-time controllability, i.e., particular gates may 
only be reachable after an extremely long time. The controllability usually fails when 



H Q is noncompact 
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34|. However, multiple control fields may greatly enhance the 
controllability, such that any gate can be achieved at arbitrary positive times (rendering 
the system strongly controllable), if their corresponding control Hamiltonians span the 
whole Lie algebra of the symplectic group. Unfortunately, many physical models possess- 
ing multiple control fields for CV quantum computing are still not strongly controllable. 
However, as discussed in Section IV, it is possible to design systems with proper coupling 
Hamiltonians that guarantee the strong controllability of the system. In what follows, 
we assume that a final time tf has been found at which the system is controllable, and 
derive the corresponding control landscape topology analytically. 

It is instructive to review the case of discrete variable quantum gate control. In this 
case, a typical measure of fidelity is the Frobenius matrix norm on U(N), defining the 
following objective function: 

J(U) = \\W - U\\ 2 F = Tr(W - U)\W - U) = 2N - 2ReTr(W^U). (11) 

where W G U(N) is the target unitary transformation. Extensive numerical studies 
of gate control optimization using this measure have been reported [3]; local iterative 
algorithms have generally succeeded in locating the global optima (albeit at a somewhat 
higher expense than the optimization of quantum observables) . Application of the tools 



of differential geometry on the Riemannian manifold U(N) 351 ] has enabled identification 
of the critical manifolds of this objective function [lOj]. The essential findings of this 
work were 1) no local traps exist in the control landscapes for discrete quantum gate 
implementation, consistent with the routine success of the associated OCT calculations 
and 2) the critical topology of these landscapes is universal, i.e., independent of the 
target gate W as well as the Hamiltonian. The correspondence of critical points on the 
unitary propagator and control field domains was also the subject of a separate recent 
investigation, where the critical topology on C(t) was derived explicitly under the electric 
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dipole approximation, and found to be essentially identical to that on the dynamical group 
U{N) (g, 36]. 



A. Critical Landscape Topology over the Symplectic Group 

The optimization problem (jSj) can be decomposed into two separate problems on the 
groups Sp(2iV, R) and M? N , respectively, where the latter has a unique optimum s = w. 
Thus, we only the critical topology on Sp(2iV, E) of the function 

J{S) =Ti{S -W) T {S -W), SeSp{2N,R). (12) 

The method to solve the critical points is to perturb the cost function in an arbitrary open 
subset and determine when the variation always vanishes irrespective to the perturbation, 
which leads to the following condition [291 ] : 



(S T S -W T S)J = J(S T S-S T W). (13) 

In contrast to the case of the unitary group, the critical topology becomes more complex 
as will be shown below, because (1121) does not possess a favorable linear form as shown in 
pip . Nevertheless, this equation is still analytically solvable owing to its highly symmetric 
form. Here we only present the conclusions, and the details can be found in a separate 



work of the authors 29|. Let W = UEV be the singular value decomposition of W, 
where U, V £ OSp(2iV, M.) and the diagonal matrix E contains the singular values of 
W. Suppose that there are 2no singular values e$ = 1 (if existing) and the remaining 
singular values 1 < e\ < ■ ■ ■ < e s have degeneracies ni, • • • ,n s . In the context of quantum 
optics, the pre-transformation V and the post-transformation U decompose the target 
symplectic transformation into operations on decoupled modes, which includes no linear 
operations and the remainder are squeezing operations with squeezing ratios ex, • • • ,e s 
on the p-components. In the following discussions, we always assume that U = V = I2N 
so that the physics can be discussed in a simple canonical coordinate system where all 
the modes are decoupled. 
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The critical submanifolds can be expressed as 291 ] : 

S* = R T DR, R e Stab(E), (14) 

where the stabilizer Stab(-E) of E in OSp(2iV, M) is defined as 

Stab(E) = {Re OSp(2N,R)\R T ER = E} = OSp(2n ) x 0{m) x ■■■0{n s ). 

The characteristic matrix D consists of different operations on the separate modes rep- 
resented by the following diagonal blocks (their inverses appear pairwisely symmetrically 
in D and E according to the reciprocal property of eigenvalues of symplectic matrices): 

(1) Type-I operations D' a = e a I m > a corresponding to a sub-block E' a = e a I m > a in E, 
which are identical with those operations in the target transformation W; 

(2) Type-II operations Dp = — e^ 3 /™" corresponding to a sub-block E'p = e^lm" in 
E, which reverse the directions of the quadrature vector of the corresponding modes and 

1/3 

is followed by a squeezing operation with ratio on the g-components. 

(3) Type-Ill operations 

jj I C<5 / &QSXygI m ^ s i sin X*y$ Irrijs 

V e 7 \ ± 

where = arccos — e i/^s- e s/e~, _ ^ ^ corresponding to a sub-block 

Ejs = . I , e/ < e 7 < ej, 

which rotates two decoupled 7-th and 5-th modes with an angle x^ followed by a uniform 
squeezing operation on both modes. 

(4) The sub-block for the particular singular value eo = 1 is 



Dr 



Im 



'I no -mo 



which leaves mo modes of linear operations invariant and reverses the direction of quadra- 
ture vectors of the other no — m modes in the phase space. 
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Moreover, the critical transformations allows for arbitrary linear optical transforma- 
tions in the stabilizer Stab(.E') as shown in (TBI) . Each combination of the above operations 
forms a critical submanifold, and the indices {m Q , m' a , m'L m^s} label each individual crit- 
ical submanifold. All admissible combinations can be enumerated to count the number 
of the critical submanifolds. 

Having identified the complete set of critical points, it is important to determine 
their optimality status (i.e, maxima, minima or saddle points) through analysis of the 
eigenvalue structure of the Hessian quadratic form, so as to understand their influence 
on the search effort for optimal controls. Enumeration of the positive, negative and zero 
eigenvalues of the Hessian at a critical point provides information about the number of the 
upward, downward and flat orientations of the principal axis directions of the landscape, 
respectively. The counting results are given in 29(; they show that the local minimum 
S* = W is unique among all the critical solutions, and the remaining critical points are all 
saddle points. Thus, no false traps are present to impede the search of optimal controls. 



B. Landscape critical topology on OSp(2A r , R) 

Carrying out optimal control field searches using only linear quantum optics corre- 
sponds to searching over only the compact subgroup OSp(2A r , R). Such a strategy may 
be desirable in the case of compact symplectic gates, such as F, X or Z above. A priori, 
it is not obvious when it is desirable to search over only this subgroup versus the whole 
group using a combination of linear optics and squeezing. Since such decisions would be 
facilitated by knowledge of the landscape topology over OSp(2iV, R), we also carry out a 
critical point analysis over this domain. 

Firstly, the landscape function is greatly simplified on OSp(2iV, R): 

J(S) = Tt(S -W) T (S -W) =4N -2Tr(W T S), S G OSp(2iV, R), (15) 

where W is the target compact symplectic gate. In fact, by (j3J), the control landscape 
over OSp(2iV, R) can be mapped to an equivalent control landscape in the form of ffTTT) . 
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Hence, it is obvious [29] that the critical topology is identical with that of the unitary 
;ransformation landscape l|. Specifically, the critical manifolds are OSp(2iV, R) orbits 
29|, i.e., S* = WO T D m O, where O G OSp(2iV,R) and 

l-I m \ 

lN~m 

3-m 



D r . 



\ 



m 



0,1, 



N. 



iN-m ) 

Hence, there are a total of N + 1 critical solutions, with values of the cost functional 



J = 0, 8, 16, 8N. The minimum and maximum values of J correspond to S = W and 
S = —W, respectively, and the rest critical points are saddle. 

We note the important fact that whereas the dimension of a unitary transformation 
representing a iV-qubit discrete quantum gate scales exponentially as 2 N jl^j], for contin- 
uous systems, the dimension of the symplectic transformation required to implement the 
equivalent continuous quantum gate within the Clifford group scales linearly as 2N, where 
iV is the number of qunits. This is understandable because the computational capabil- 
ities of discrete quantum computers extend beyond those of Clifford-gate CV quantum 
computers. 



IV. OPTIMAL CONTROL OF CONTINUOUS VARIABLE GATES 

Algorithms for quantum optimal control have been extensively developed for the con- 
trol of discrete (finite-dimensional) quantum systems, or continuous quantum systems 
that can be treated as finite-dimensional to a reasonable approximation. Broadly speak- 
ing, two distinct types of discrete variable quantum control problems have been consid- 
ered from an algorithmic point of view: 1) control of quantum observables; 2) control 
of dynamical propagators (gates). It has been noted that the latter is generally compu- 
tationally more expensive, in part because the solution set is relatively small to locate. 
One of our primary aims in this work is to characterize the expense of solving 3) optimal 
control problems for CV quantum dynamical propagators under various constraints and 
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to compare this expense to that of 2). For similar reasons, one would expect problem 3) 
to be inherently more difficult than the optimization of CV observables. 

The analysis above shows that despite the noncompactness of the symplectic group 
of CV propagators, it is possible to choose a scalar fidelity function with a fairly simple 
critical topology. In this section, we describe the physical models employed for CV 
quantum computations and the numerical algorithms used to search for optimal controls. 



37]. 



A. Physical control systems for CV quantum information processing 

Over the past few years, several physical models have been suggested for the im- 
plementation of symplectic gates. Early proposals focused on coupled pairs of co njugat e 
continuous variables describing quadrature modes of the electromagnetic field lj], Q, 
In such traditional quantum optics models, it is typically possible to apply only one con- 
trol at any given time. Although it is often possible to implement CV gates via one control 
with reasonable fidelity if the final time tf is chosen judiciously, the quality of control will 
be downgraded. In the present work, we do not limit ourselves to these restricted control 
Hamiltonians from quantum optics, since our goal is to explore the generic properties of 
the CV gate optimal control problem. 

More recently, CV models displaying greater flexibility have been proposed that may be 
more suitable for optimal control of quantum gates. In particular, these control systems 
allow for the simultaneous application of two or more independent controls. A simple 



model raised in 
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38] considered the off-resonant interaction of light with a collective 



spin described by the effective Hamiltonian 

which represents a strongly coherent light beam polarized along the x-axis that propagates 
along the z-axis through the atomic ensemble. In addition, it is assumed that arbitrarily 
fast local phase shifts are implementable by single-model control Hamiltonians: 

Hi = x\ + pi, H 2 = x\+ p\. 
16 



This system satisfies the rank condition, but is not ensured to be fully controllable because 
Hq is noncompact. In 



371 ]. control pulses are restricted to be instantaneous and exerted 
in certain sequences to simplify the analysis and the experimental realization. However, 
here we will assume that the control pulses can be arbitrarily shaped, so that a greater 
degree of precision is possible in tailoring the control Hamiltonian to match theoretical 
predictions. 

Recently, atomic ensembles, particu 
a promising medium for CV QIP 
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arly ensembles of trapped ions, have emerged as 



40, 41] , because the trapped ions are thermally 



isolated from their environment, minimizing decoherence effects. These systems may also 
offer a degree of flexibility suitable for gate synthesis via OCT. For CV gates, quantum 
information is stored in the vibrational modes of the trapped ions. To couple (entangle) 
these vibrational modes, several studies have examined the interactions of vibrational 
states of trapped ions with some quantized fields inside an optical cavity 0, 4ol |. through 
which it is possible to indirectly tune the coupling between vibrational modes. 

For concreteness, consider a model wherein two trapped ions with internal electronic 
levels are coupled to external lasers. They are also coupled to a cavity mode with fre- 
quency u c , described by annihilation and creation operators a and a\ respectively, where 
the harmonic frequency of each trap is 2v. We assume that the cavity is oriented along 
the x axis and the laser is incident along the y or z axis. Then in a frame that is rotating 
with frequency u c , the interaction Hamiltonian coupling the vibrational modes of the ions 
with the cavity and laser fields can be written (omitting the electronic states): 

H = 2v{b\b x + b\b 2 ) + V, 

where bj and fej (j = 1, 2) are the annihilation and creation operators of the vibrational 
modes. The interaction Hamiltonian V is a function of the coupling constants between 
ions and lasers and single-photon coupling strengths. Under reasonable assumptions 
regarding the size of the traps compared to the laser wavelength, and with proper detuning 



of the lasers from the cavity mode, it can be shown 40j that the interaction Hamiltonian 
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for the first ion can be approximated as 

V w ru^&i + ah\) + r 2 i{(Jb\ + ah)- 

where the r mn are functions of the frequencies and coupling constants. This system 
actually involves three harmonic oscillators, where ion-cavity interaction Hamiltonians 
produce indirect couplings between the vibrational modes of the two ions. By modulat- 
ing the frequencies of these lasers through time, we can achieve time-dependent control 
Hamiltonians necessary for the implementation of CV gates with optimal fidelity. The 
controls in this model induce nonlocal interactions between qunits. The associated control 
system does not satisfy the controllability rank condition, and hence is uncontrollable. 

The above models display features that are representative of current proposals for CV 
QIP. As we will see below, their controllability properties are of particular importance. 
The light-collective spin interaction control system (hereafter referred to as the "photon 
model") is sufficient for achieving arbitrary symplectic transformations in experiments, 
but its controllability is not strong enough to achieve the target in arbitrary finite time. 
On the other hand, since the ion trap model is not fully controllable, for certain gates, 
there does not exist a final time tf at which the gate can be reached with arbitrary 
precision. Nonetheless, for many gates, it is possible that such a final time exists. 

In addition to these two physical models, we also choose a strongly controllable system 
below in order to explore the effects of controllability on the properties of gate optimiza- 
tion. The following control Hamiltonians were employed for this purpose: 

H\ = a{ 2 — a\ + a} 2 2 — a\, H 2 = d\di — a\a2 + i(a\ + di){a\ — a 2 ); 

the internal Hamiltonian in this case consisted of uncoupled harmonic oscillators, i.e. 
H = a[di + d\a2 + 1. 

As in the case of discrete variable quantum control, it is possible to impose additional 
constraints on the optimization problem based on the physical implementation of choice, 
such as bounds on the control intensity or the time derivative of the control pulse. A 
study of the impact of such constraints is properly the subject of a separate work. 
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In the following, we make several comparisons to discrete QIP. For this purpose, we 

n 

assume the standard physical model of NMR-based quantum computation [6j. In this 
model, the internal Hamiltonian H consists of nuclear spins that are coupled in the 
absence of the control field. The coupling between N spins (only up to 2-qubit interactions 
are considered) is achieved through standard NMR coupling Hamiltonians of the form: 

N N N 

i=l i ,j=l a,/3=x,y,z i=l 

where a^' y,z are the standard Pauli matrices representing observables of the i-th qubit. 
The first term splits the energy levels via a static magnetic field along the z-axis, with 
uji being the Raman frequencies; the second term represents the internal couplings be- 
tween the qubits (e.g., chemical shifts); the last term, the control Hamiltonian, interacts 
each qubit states to a time- variant x-axis radio-frequency control field Ci(t). Because 
of the tensor product structure of qubit subspaces in these expressions, the total system 
dimension scales as 2^. One can verify that this system is controllable, but not strongly 
controllable. 



B. Numerical implementation 

Several optimization algorithms, including iterative methods such as the Krotov algo- 



rithm [3|and tracking methods such as D-MORPH 42j, have been applied in the OCT 
of discrete quantum gate implementations. These algorithms can vary considerably in 
optimization efficiency, but they are all based on information pertaining to the first func- 
tional derivative of the objective function with respect to the control field. Since our 
primary goal in this work is to compare the properties of discrete and continuous gate 
OCT, we adopt gradient algorithms to search for optimal controls, which, although not 
the most efficient, offer the simplest and most direct opportunity for comparison. In 
particular, they are ideal for detecting landscape traps. By comparing the magnitudes 
of the gradient and Hessian of the objective function along the search trajectory, we can 
obtain a understanding of the factors that govern optimization efficiency and what can 
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be applied to the design of tailored algorithms in future work. 

The electric field Cj(s, t) was stored as a p x q matrix, where p and q are the number 
of discretization steps of the algorithmic time parameter s and the dynamical time t, 
respectively. For each algorithmic step Sk, the field was represented as a q- vector for the 
purpose of computations. Starting from an initial guess Ci(so,t) for the control field, 
the equations of the motion were integrated over the interval [0,tf] by propagating the 
Schrodinger equation over each time step i& — > tk+i, producing the local propagator 
U(tj+i,tj) = exp[—iH(si,tj)tf/(q — 1)]. The method used for this purpose differed for 
the discrete quantum and symplectic gate optimization problems. For discrete quan- 
tum systems, the Hamiltonian matrix was diagonalized (at a cost of iV 3 ), followed by 
exponentiation of the eigenvalues, and multiplication of the resulting matrix on the left 
and right by the matrix of eigenvectors and its transpose. Alternatively, a fourth-order 
Runge-Kutta integrator can be employed for the propagation. The propagation toolkit, 
which involves precalculating the matrix exponential at discrete intervals over a specified 
range of field amplitudes, was used to further improve the speed of Hamiltonian integra- 
tion for discrete quantum systems. The initial guess C(so,t) was alternatively taken as 
a random, constant, or sin pulse field. 

The Pade approximation for the exponential function was used to calculate the local 
symplectic propagators S(tj +1 , tj) = exp[Jif (sj, tj)tf/ (q — 1)]. Since the matrix JH(si, tj) 
is not symmetric, it is not possible to calculate its exponential via diagonalization and 
subsequent scalar exponentiation of its eigenvalues. The type (p, q) Pade approximation 
for e x is the (p, g)-degree rational function P pq (x) = N pq (x) / D pq (x) obtained by solving 
the algebraic equation YlkLo xk ~ N pq (x) / D pq (x) = 0(x p+q+1 ), i.e., P pq (x) must match 
the Taylor series expansion up to order p + q. The primary drawback of the Pade approx- 
imation is that it is only accurate near the origin, so that the approximation is not valid 
when JH(si,tj) is too large and when its eigenvalues are not too widely spread. For the 
problems considered in this work, the norm of Jif (sj, tj) was always small enough for high 
accuracy in the approximation. Because of the noncompactness of the symplectic group, 
which results in the functional derivatives of the objective function changing too rapidly, 
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Runge-Kutta integration was not used. Due to the large number of iterations generally 
required for convergence of CV quantum controls, the speed of the matrix exponentia- 
tion algorithm is particularly important. However, the propagation toolkit was generally 
found to be inadequate for speeding up symplectic matrix propagation; discretization of 
the control field amplitudes produced unacceptable errors in the matrix exponential, and 
the maximum amplitudes often grew abruptly during optimization. 

Minimizations of the fidelity function were typically done using the Polak-Ribere vari- 
ant of the conjugate gradient (CG) method. Step size was varied adaptively based on 
Brent's method for line minimizations. In several cases, adaptive step size steepest de- 
scent was employed to analyze the behavior of gradient flow lines. Steepest descent was 
found to converge much less efficiently than CG for the symplectic gates. Its performance 
for discrete unitary gates was considerably better. For these algorithms, the gradients of 
the objective functions were calculated analytically via the following expressions: 



where H,i(t) = W{t)HjU{t), for discrete unitary gates. Here, the H^s are the Hamiltoni- 
ans that couples to the time-dependent control; for molecular systems, it represents the 
dipole moment operator. 

Integration of the gradient flow equations on the domain of symplectic or unitary 
propagators was carried out using the fifth-order adaptive step size Runge-Kutta method. 

C. Impact of landscape topology on numerical control search 

From the critical topology analysis, it is clear that optimal control landscapes for sym- 
plectic transformations are devoid of suboptimal local traps, regardless of the structure 
of the gate transformation or the Hamiltonian of the system. Additionally, all critical 




Tr [H i (t)J(S T (t f )W-S T (t f )S)} 
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points stay in a bounded region in the symplectic group, i.e., they can affect the optimal 
search only when it enters that region. 

The compactness and degeneracies of the singular values of W determine the critical 
topology of the control landscape. When W is compact (e.g., F, X and Z or combinations 
thereof), the number of critical manifolds in the control landscape scales linearly as N + l, 
where N is the number of qunits. When W is noncompact and has fully degenerate 
singular values (i.e., = e* ^ 1), the number of critical submanifolds can be shown to 



be quadratic in iV 29]: 



N= (iV + 2)V2, iVeven; 
\ (iV+l)(iV + 3)/2, N odd; 

The scaling for nondegenerate gates shoots up when the degeneracy is broken. For the 
case that W has fully non-degenerate singular values, the upper bound for the number 
of critical submanifolds is 

'" 2 N - 3m N\ 



■f m\(N-2m)V 
which is super-exponential. 

By contrast, it was previously shown that the number and critical values of all discrete 
quantum logic gates are independent of the gate and depend only on the dimension of the 
system jjj. FigJT] compares the scalings calculated above with that of (iV-qubit) unitary 
gates. In most cases, the number of critical submanifolds for symplectic gates grows 
much faster than that for unitary gates. However, the number of critical submanifolds is 
always finite, and hence the critical region is bounded for arbitrary target gates W, and 
contained in the ball centered at W with radius 



R 



N 



v 5> 2 + e- 2 + 3ef + 3e7^) > 2V2N 



equal to the distance to the farthest critical points. The volume of this region is roughly 
of the order V ~ R 2N2+N . Assume that the attraction is effective in a r-ball around each 
critical point; then the ratio o of the volume of the attractive regions to the volume of 
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the region of critical points 

^f xr 2N 2 +N , r ^2N 2 +N 

can be numerically proven to go rapidly to zero when N approaches to infinity for ar- 
bitrary r < R, impling that the probability of a random initial guess starting close to a 
saddle manifold should be so small as to be negligible in practice. 
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Fig 1: The scaling of the numbers of critical submanifolds for fully non-degenerate (iV-qunit) 
symplectic gates and (iV-qubit) unitary gates. 

For a compact target symplectic gate, the critical submanifolds are completely iden- 
tical when the search is carried out on the full symplectic group (corresponding to linear 
optics plus squeezing) or its compact subgroup (corresponding to only linear optics), the 
only difference being an increase in the dimension of the search space. Therefore, the 
additional directions accessible through squeezing transformations are not expected to 
improve convergence toward the optimal solution. These observations collectively paint 
a fairly simple picture of CV gate landscape topology, which, although more compli- 
cated than the topology of discrete gate landscapes, should not preclude efficient control 
optimization. 

Throughout this paper, we will use gradient algorithms to optimize the control field. 
The search process in the kinematic picture can be represented by the so-called gradient 
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flow on the symplectic group: 



dS 



-VJ(S). 



It is instructive to estimate the convergence speed of the gradient flow, via linearization 
of this equation in a sufficiently small neighborhood of the global optimum S = W, which 
gives 



d SS 



-Hess(jy) SS, 



where SS is the deviation of S from W, which is proportional to the gradient vector at 
S. The positive definite matrix Hess(U/) is the Hessian matrix at W. Therefore, the 
convergence of SS is exponential and its rate is dominated by the smallest eigenvalue of 
Hess(U/), identified as e~ 2 < 1. By comparison, a similar estimate for gate search on 
the discrete unitary group reveals a constant convergence rate of 1. Therefore, search for 
noncompact symplectic gates will display slower convergence in general, depending on 
the magnitudes of the singular values of the target gate. In particular, the convergence 
speed for the phase gate (or squeezing gate) decreases with increasing phase shift (or 
squeezing ratio). Fig J2] shows the convergence speed of gradient flows for the SUM gate 
on the symplectic group, and, for comparison, the CNOT gate on the unitary group. 
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Fig 2: The convergence of the gradient flows for optimal search of the SUM gate on Sp(4, M) 
and CNOT gate on W(4). 
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V. NUMERICAL SIMULATIONS 



In this section, we numerically solve for optimal control fields for Clifford group gates 
and composite CV algorithms, and compare the optimization effort and control field 
complexity to those of the corresponding gates over discrete variables. As such, we aim 
to identify distinguishing features of CV gate control that will have the greatest impact 
on its computational and experimental implementation. 

A. SUM Gate 

Consider the control problem where W is the SUM operation, whose matrix form is 
shown in ([7]). We simulate the realization of the SUM gate using both the photon model 
(k = 1) and the strongly controllable system described above. Figj3] compares the effects 
of strong versus weak controllability on convergence speed, starting from either a random 
guess or near saddle solutions. Although the SUM gate is reachable with both models at 
the final time, the weaker controllability of the photon model compromises convergence 
speed. In addition, FigJH shows that the corresponding optimal control fields are more 
expensive in that their fluences are much greater. 

The singular values of SUM are E = diag{l. 618, 1.618, 0.618, 0.618}. The analysis 
above predicts that there should be 4 critical submanifolds for a degenerate 2-qunit 
gate. The first one D\ = diag{l. 618, 1.618, 0.618, 0.618} contains one type-I block, 
whose corresponding critical submanifold is the isolated global minimum SI = SUM. 
The second, D2 = diag{— 0.852, —0.852, —1.174, —1.174}, contains one type-II block, 
whose corresponding critical submanifold is an isolated saddle point. The third, D 3 = 
diag{l. 618, —0.852, —1.174, 0.618} contains one type-I and one type-II block. The last 
saddle contains a type-Ill block, which actually contains two isolated points. In sum- 
mary, there are a total of 5 critical submanifolds, including 4 isolated points and 1 one- 
dimensional manifold. The Hessian analysis is summarized in Table I. 
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No. 


Critical value 


D 


D+ 


D- 


type 


1 








14 





minimum 


2 


18.623 





10 


4 


saddle 


3 


9.311 


1 


12 


1 


saddle 


4(±) 


10 





11 


3 


saddle 



Table I: Critical topology for the SUM gate. 



As discussed above, the saddles will rarely be encountered during the progress of most 
optimization trajectories. This is also supported from the simulation result (Figj3]), the 
saddle manifolds appear to have a slightly adverse effect on optimization efficiency for 
the strongly controllable system, and almost have no influence on the convergence of the 
photon model. 

Note that the free Hamiltonian for the photon model happens to be proportional to 
the matrix logarithm of the SUM gate; thus, using this model, SUM can be achieved 
merely via free evolution in k -1 units of time. Simulations show that the SUM gate is 
always realizable in any time longer than k~ 1 (e.g. FigHJ). However, it is interesting to 
see if the local controls permit the gate to be achieved in a shorter time. Figj5] shows an 
example employing a final time less than k~ 1 . As can be seen, the control search does not 
converge. By contrast, for the strongly controllable system, optimal control fields exist 
even for very small tf, although the expense increases and the shape of control fields 
tends to become more singular (Fig|6j) . 

An important feature not shown in the figures is that the optimal search suffers from 
serious numerical instabilities when it starts far away from W because the system dynam- 
ics involves exponentially increasing components. This can easily occur when there is no 
a priori knowledge of an appropriate choice for the control fields. By contrast, because of 
the compactness of the dynamical group, the control optimization for discrete quantum 
systems never encounters this problem. 
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Fig 3: The convergence of dynamical search for the SUM gate using conjugate gradient algo- 
rithms with the photon model (left) and a strongly controllable system(right). 

B. SWAP gate 

Optimal control landscapes over 0Sp(2iV, R) and U(N) have identical topology and 
geometry, because these dynamical groups are isomorphic. This section simulates the 
optimal search of a (compact) transform that swaps the states of the two qunits as 
follows 

/o 1 o\ 

w= 1000 . 

1 

yO 1 Oy 

For this gate, we employ the photon model with two local phase controls, using two 
different free Hamiltonians: (i) H = X1P2, which involves squeezing operations; (ii) 
H = X1P2 — X2P1, which involves only linear optics. It can be verified that these systems 
are controllable over (i) Sp(4, R) and (ii) OSp(4, R), respectively. The simulation results 
in Fig IT] show that the optimal search restricted on OSp(2iV, R) generally exhibits fast 
convergence, as in the case of control of discrete unitary gates. This is not surprising 
because the group OSp(2iV, R) is isomorphic to the unitary group 14 (N) and the corre- 
sponding dynamical control system is equivalent to a iV-level discrete quantum control 
system. By contrast, optimal control using squeezing operators as well exhibits no ad- 
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Fig 4: The optimal control fields for CV SUM gate control after searching from a random initial 
guess, for (a) the photon model and (b) a strongly controllable model. 

vantages compared to using only linear operations; in fact, the resulting control fields 
have much greater ffuences (FigJHJ). 

Because OSp(2iV, R) is isomorphic to U(N), a comparison of the gradients of the fi- 
delity on Sp(2iV, R) and OSp(2iV, R) sheds light not only on the comparative efficiencies 
of these two optimization problems, but also the origin of the slower convergence of non- 
compact CV gate optimization versus that of discrete gate optimization. Figj9] displays 
the norm of the gradient of the fidelity of the SWAP gate with respect to the control 
field on Sp(2iV, R) and OSp(2iV, R), at each algorithmic step during the course of opti- 
mization. In order to sample more points along the optimization trajectory, a steepest 
descent algorithm was employed in this case, starting from near a saddle point of the con- 
trol landscape. As can be seen, the norm of the gradient is larger on OSp(2iV, R) at most 
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Fig 5: Optimal control search for the SUM gate using conjugate gradient algorithms with the 
photon model at a final time smaller than k^ 1 . 
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Fig 6: Dynamical search for the SUM gate using conjugate gradient algorithms with the strongly- 
controllable system model at a small final time. 

points along the optimization trajectory. In addition, over several runs, it was found that 
the gradient norm changes more abruptly during optimization on Sp(2iV, R) compared 
to 0Sp(2iV, R) (or, equivalently, U(N)). It was also observed that the components of 
the gradient at successive dynamical time points change abruptly on Sp(2iV, R). These 
features, presumably originating in the noncompactness of Sp(27V, R), undoubtedly act 
to retard the convergence of searches carried out on this group. 
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Fig 7: The convergence of optimal search for the SWAP gate with linear and squeezing couplings 
(Sp(4,R)) and linear couplings (OSp(4,R)). 

C. Composite algorithms 

Composite algorithms composed of large numbers of low-dimensional Clifford group 
gates will generally have richer structures in their singular values and a greater number 
of critical manifolds. The geometry of the control landscape is also expected to be 
complexified for such gates. In practice, representing composite Clifford group algorithms 
through sequences of 1-qunit or 2-qunit gates is generally preferred. However, since using 
a larger number of gates may increase the likelihood of information loss through quantum 
decoherence, the implementation of higher dimensional transformations is desirable in 
some instances. 

Consider a composite operation on 3 qunits that sums the values of their g-components. 
This gate can be decomposed into two elementary gates SUM(1, 2, 3) = SUM(2,3) x 
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Fig 8: The control fields for the swap gate with (a) linear and squeezing couplings (Sp(4, ] 
and (b) linear couplings (OSp(4, K)) . The solid lines refer to the control field on the first qunit, 
while the dotted refer to the second. 
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This composite operation can be implemented using a photon model where the internal 
Hamiltonian includes interactions between 1-2 and 2-3 qunits, i.e., Hq = X1P2 + X2P3, and 
three local control Hamiltonians applied in the form of Hj = x 2 - + p?, j = 1, 2, 3. Again, 
because the internal Hamiltonian H is noncompact, full controllability is not guaranteed 
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Fig 9: The norm of the gradient vector during optimal search for the SWAP gate, starting 
from near a saddle point, with linear and squeezing couplings (Sp(4, R)) and linear couplings 
(OSp(4, R)). The left panel is for the first control field and the right panel is for the second. 

at arbitrary final time tf. In the case of the ion trap model, the internal Hamiltonian H 
consists of uncoupled harmonic oscillators, and two nonlocal controls 

H\ = a)bi + ab\, Hi = a)b 2 + abi 

are applied. As discussed above, there is no guarantee that the gate will be reachable at 
any final time within this model, since the system does not satisfy the controllability rank 
condition. In FigflQi we compare the convergence speeds of optimal control search for 
the 3-qunit SUM gate using these models with that of its discrete quantum counterpart, 
the Controlled-CNOT gate (Toffoli gate). In the latter case, the NMR control system 
described above was used. 

As can be seen from Fig. 10, neither the ion trap nor the photon control search converges 
within the specified tolerance for the chosen final time, whereas the strongly controllable 
system does converge. The weakly controllable and uncontrollable systems therefore 
display similar behavior for this composite gate. This example demonstrates that the 
controllability of a CV gate control system may become a more important consideration 
for higher dimensional gates. 

For both discrete and continuous quantum systems, the decrease in convergence speed 
with increasing system dimension appears to be severe; in addition, control field searches 
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Fig 10: The convergence of optimal searches for the 3-qunit SUM gate with photon model and 
ion-trap model, and 3-qubit Controlled-CNOT gate with NMR model. 

are more likely to become trapped due to easier loss of controllability. Therefore, for a 
quantum algorithm involving a polynomially large number of operations (i.e., primitive 
gates), it would indeed appear more efficient to apply sequences of smaller gates rather 
than attempting global search over transformations on the whole set of qunits. The 
scaling of optimal control search effort with system size for CV gates is an important 
subject for future study. 

Finally, from FigJTT} we observe a stark difference in the fluence of the optimal control 
fields obtained using local versus nonlocal controls. For the photon model (local controls), 
the fluence of the optimal fields exceeds physical limits, whereas for the ion trap model 
(nonlocal controls) and the NMR model, the fluence remains bounded. This indicates that 
for composite gates, CV control models employing nonlocal controls may be preferable 
to those whose design requires the use of local controls. Note that although the NMR 
model employs local controls, their fluence remains small, suggesting that this problem 
does not arise for discrete variable gates. 



33 




Fig 11: The optimal control fields for (a) 3-qunit SUM gate with ion trap model, (b) 3-qunit 
SUM gate with photon model and (c) 3-qubit Controlled-CNOT gate. For (b) and (c), the solid 
(dashed, dotted) lines are the control fields on the first (second, third) qunit/qubit. 

VI. DISCUSSION 

The absence of local traps in control landscapes for symplectic gate fidelity indicates 
that given sufficient time, local gradient-based algorithms will generally succeed at reach- 
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ing the global optima (perfect fidelity), assuming the system is controllable. This prop- 
erty, combined with other attractive features of continuous QIP such as the high band- 
width of continuous degrees of freedom, strengthens the feasibility of QIP over continuous 
variables. 

We have seen that CV gate optimization problems can be divided into two classes 
with inherently different complexities. The first, wherein only linear Hamiltonians are 
employed as controls, is mathematically identical to the problem of discrete unitary gate 
optimization. The second, in which squeezing operations are also employed, is generally 
more expensive. 

A second point of distinction between these two control problems is the complexity of 
the optimal fields. The optimal fields for discrete gate control are typically in resonance 
with the transition frequencies of the system, within the weak-field regime. By contrast, 
in the case of CV gate optimization over Sp(2iV, R), the fields are usually not simply re- 
lated to the natural resonant frequencies of the system because the Hamiltonian possesses 
imaginary eigenvalues. These eigenvalues produce exponentially increasing and decreas- 
ing modes in the field. The former can result in instabilities during the optimization 
process 44j]. Moreover, for CV gate optimization over Sp(2iV, K), there is a strong depen- 
dence of field complexity on the final dynamical time. It is important to identify several 
controllable tf's and choose the one corresponding to control fields that are physically the 
most simple to implement. However, the higher degeneracy of these fields means that we 
are presented with more choices for convenient physical implementation, and point to a 
rich variety of distinct control mechanisms that reach the same objective. 

Finally, the effects of quantum system controllability on control optimization are more 
subtle for CV gate control than for discrete variable gate control. For CV gate control, 
it is often difficult to identify a final dynamical time tf at which the gate can be reached 
with high fidelity, if the system is not strongly controllable. Moreover, several current 
physical models for CV gate synthesis are not fully controllable for any choice of tf. For 
these systems, even when the gate of interest is reachable, the cost of optimal search is 
typically steep. It is therefore particularly important to employ physical control systems 
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that satisfy the conditions for controllability. 

These details must be borne in mind when assessing the comparative difficulty of im- 
plementing quantum communication protocols over continuous versus discrete variables, 
and underscore the importance of using optimal control theory in the design of high- 
fidelity symplectic gates. An important conclusion of this work is that it is important 
to use the methodology of control theory when designing physical systems for the im- 
plementation of continuous variable quantum gates, both in the choice of appropriate 
physical systems and in the determination of the controls themselves. Application of 
OCT to practical CV quantum information tasks will require the imposition of penalty 
terms on the control Hamiltonians corresponding to physical constraints. For control of 
discrete states in molecular systems, the most significant constraint on the control fields 
is the total fluence, since current pulse shaping technologies are capable of producing the 
majority of shapes predicted by OCT. Ongoing efforts to design pulse shapers with ultra 
high bandwidth should further facilitate implementation of the theoretically predicted 
fields. By contrast, for CV systems, it is difficult to shape control pulses within certain 
physical models. Applying shape constraints to the optimization problems above would 
amount to choosing amongst the highly degenerate sets of control Hamiltonians that solve 
the generic OCT problem. 

In the context of particular Hamiltonians, time optimal control theory has been applied 
to assess the minimal time necessary to implement various primitive discrete quantum 
gates A natural counterpart to the current work is the application of time 

optimal control theory to symplectic gates, using restricted control Hamiltonians that 
are commonly implemented in the quantum optics laboratory. Given the susceptibility of 
CV quantum systems to noise, time optimal control of CV gates is particularly relevant. 
Such studies would necessitate determination of the minimal length geodesies in the 
noncompact symplectic group, or in subgroups of Sp(2iV, K) that can be reached using 
the restricted controls. Preliminary work [26] along these lines employing the so-called 
Bloch-Messiah theorem has been reported, but it is unlikely that analytical solutions 
exist for most time optimal symplectic transformation control problems, especially for 
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higher qunit systems. As such, the development of numerical OCT algorithms suited to 
time-optimal CV gate control is an important future challenge. 

In the present work, we have adopted the approach of optimizing a scalar objective 
function for gate fidelity. An alternative approach would be to track a predefined trajec- 
tory on the space of symplectic propagators between the initial and target gate. From 
a local perspective, the former is computationally more efficient, but globally, the latter 
may offer an advantage. Perhaps more importantly, the latter approach would provide 
insight into the effect of the map between control fields and dynamical propagators on 
optimization efficiency. Indeed, the numerical simulations above indicate that the prop- 
erties of this map are primarily responsible for the differences in optimization efficiency 
between discrete and CV quantum OCT. Moreover, these properties will impact the ef- 
ficiency of control optimization for any observable of a quadratic CV system, since the 
symplectic transformation uniquely defines the infinite-dimensional unitary propagator 
corresponding to the CV quantum dynamical process. 

Finally, we note that universal quantum computation requires nonlinear symplectic 
gates corresponding to the ability to count photons in the electromagnetic field [43]. The 
implementation of nonlinear symplectic gates necessary to achieve universal continuous 



quantum computation is known to be difficult to achieve with high fidelity [25|]. Purifica- 
tion protocols are necessary to distill from an initial supply of noisy nonlinear symplectic 
states a smaller number of such states with higher fidelity. Future work should consider 
the challenges inherent in implementing such gates through the methodology of optimal 
control theory. 
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